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I. INTRODUCTION 

Higher-dimensional gravity appears in the landscape of 
modern physics as a result of the string theory or in var- 
ious brane-world models. Both, string theory and brane- 
world contexts imply the existence of higher dimensional 
black holes. For the four dimensional gravity, there exists 
a uniqueness theorem that states that the asymptotically 
flat stationary black hole is described by the well-known 
Kerr solution. In higher dimensions, such a uniqueness 
theorem does not hold any more, but a great variety of 
solutions describing black objects, such as black holes, 
strings, branes, rings, saturns, exist instead There- 
fore, the stability of rotating black holes in higher than 
four space-time dimensions is an open appealing prob- 
lem 0: 0- Stability could help to select among the va- 
riety of "black" solutions those which could be realized 
in nature. Another motivation to study the evolution of 
higher-dimensional black hole perturbations comes from 
the possibility to observe mini-black holes in future ex- 
periments with cosmic showers or in the Large Hadron 
Colhder [2|. 

The stability of spherically symmetric static black 
holes has been intensively studied in recent few years. 
It started from the papers by Ishibashi and Kodama 
El I 0i 0i who, using the gauge- invariant formalism 
5| , proved that the D-dimensional generalization of the 
Schwarzschild black hole, described by the Tangher- 
lini metric, is stable. Then, the stability of higher- 
dimensional black holes allowing for a cosmological con- 
stant was shown in 9] for D = 5, 6, ... 11 Schwarzschild- 
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de Sitter (SdS) and in [T^ for Z) = 5, 6, ... 11 Reissner- 
Nordstrom-anti-de Sitter (RNAdS) black holes. The D- 
dimensional Reissner-Nordstrom black hole was shown to 
be stable, although Reissner-Nordstrom-de Sitter black 
holes proved to be unstable in higher than six space-time 
dimensions [11], when both the charge and A- term are 
large enough. The stability of black holes in higher di- 
mensional generalizations of the Einstein gravity, such 
as the Gauss-Bonnet and Lovelock theories, was consid- 
ered in [l3], dill H, and an instability for large val- 
ues of the Gauss-Bonnet coupling a was found. This 
Gauss-Bonnet instability was shown to have quite pecu- 
liar behavior: the instability starts dominating after a 
long period of damped quasinormal oscillations. In addi- 
tion, the stability analysis and time domain evolution of 
perturbations were considered for the Kaluza-Klein type 
solutions, such as (squashed) Kaluza-Klein black holes 
or black strings [l6j . 

While the stability of static black holes in higher di- 
mensions is relatively well studied, the stability of ro- 
tating black holes for D > 4 is still an open question, 
because the separation of variables for all types of per- 
turbations is a complicated problem. Nevertheless, some 
particular results have been obtained for the Myers- Perry 
black holes with all equal angular momenta, though only 
for either tensor-troe perturbations pJi] or for the so- 
called zero mode [19j. 

In our previous paper f^ll] , we have started the investi- 
gation of gravitational perturbations of a much more in- 
teresting type of higher-dimensional rotating black holes, 
simply rotating black holes. These black holes have only 
one non-vanishing angular momentum component and 
are rotating in a single two-plane. Thus, they can be 
interpreted in brane-world scenarios as black holes cre- 
ated in particle collisions occurred on our brane. The 
perturbation equations for the D > 6 dimensional sim- 
ply rotating black holes, allowing for a cosmological con- 
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stant, were reduced in [2l| to the wave-like form, and 
gravitational spectrum of asymptotically- AdS solutions 
was considered. The asymptotically AdS rotating black 
holes proved to be gravitationally unstable on superradi- 
ant modes only, and with a tiny instability growth rate 

In the present paper we shall continue research that 
started in [2l| and consider tensor-type gravitational 
spectrum of asymptotically flat, Myers-Perry black holes. 
The quasinormal modes beside their importance for sta- 
bility proof, may have useful implications for the study 
of the evolution of mini higher-dimensional black holes, 



because quasinormal modes dominate at late time of the 
holding phase ^] of an evaporating black hole. 

Using the Frobenius method for the angular part of the 
perturbation equations and both WKB and Frobenius 
methods for the radial part, we have found the quasinor- 
mal modes of simply-rotating Myers-Perry black holes. 
The paper is organized as follows: Sec. II gives basic 
information on simply-rotating MP black holes. Sec. Ill 
is devoted to description of the numerical methods that 
we used for finding the quasinormal modes. Sec. IV dis- 
cusses the obtained numerical results. 
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II. BASIC FORMULAS 



The metric for the higher-dimensional Myers-Perry black hole can be written in the Boyer-Lindquist coordinates as 
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2M 
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where n = D — A, dfl'^ is the metric of the 7i-dimensional unit sphere, and 



A := (r2 + a^) 



2M 



p2 •= _|_ pQg2 , 



(2) 



Here M is the mass of the black hole and a is the angular momentum per unit mass. 

The equations for the tensor-type gravitational perturbations of this black hole are reduced to the coupled equations 
[m for the angular part 



2(1 - x^)S"{x) + (-1 + n - (3 + n)x)S'{x) +{^ + 
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x = cos{29), n — 3,4,5. 



2,3,4..., m = 0,±l,±2.... 



and for the radial part 
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with respect to the frequency uj and the separation constant Here, we have defined the tortoise coordinate r^, as 

+ 
A r) 
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III. NUMERICAL METHOD FOR FINDING OF in [111. Namely, the second-order differential equation 
THE QUASINORMAL SPECTRUM ^ is reduced to an algebraic equation with an infinite 

A. The angular part 



The differential equation for the angular part ^ can 
be solved numerically in the same way as it was done 
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continued fraction [2l|. This equation can be solved nu- 
merically to determine the separation constant fi for each 
value of uj. Thus, we can find numerically the function 
fij(io), where = 0, 1, 2 ... is an integer number. 

The number j enumerates the eigenvalues of ^. If the 
black hole does not rotate (a — 0), the equation for the 
angular part can be solved exactly to yield [2l| 

(2j + £ + \m\){£ + 2j + \m\+n+ 1), (6) 

which coincides with the eigenvalue for the harmonic 
function on the unit {n + 2)-sphere, with the multi-pole 
number (2j +£ + \m\) = 2, 3, 4 . . .. 

If the rotation parameter a is non- vanishing, the eigen- 
values II are non-integer and complex. For any value of 
w, we can enumerate them by the non- negative integer j 
in the following way. 

1. We start from the non-rotating black hole and find 
exactly the value of fi for the corresponding j . 

2. We increase the rotation parameter by a very small 
value and search for the closest to the previously 
found solution for /i. 

3. We repeat the previous step until any required 
value of a is reached. 

In this way we are able to enumerate the values of the 
separation constant for any uj and a. In order to de- 
crease the computation time, in practice, when searching 
for quasi- normal modes, we start from the known values 



of the quasi-normal frequencies for a non-rotating black 
hole [27| and, increasing step by step the value of a, find 
the closest pair of fj, and ut to the previously found one. 
An alternative approach for the computation of the an- 
gijar separation constant has been recently suggested in 



B. The radial part 

In order to simplify our notations, we will parameterise 
the black hole mass by the black hole horizon radius r+ : 

2M = rl-\rl+a^). (7) 

In solving the radial part of the perturbation equation, we 
use two alternative approaches: the numerical Frobenius 
method [23| and the semi-analytical JWKB approxima- 
tion [Hi. 

Let us start from the description of the Frobenius 
method. We are searching for the solution of the equa- 
tion Q with the quasi-normal boundary conditions that 
requires pure in-coming waves at the event horizon and 
pure out-going waves at spatial infinity. For a black 
hole in the asymptotically fiat background, they imply 
a purely ingoing wave at the event horizon and purely 
outgoing wave at the spatial infinity. If P{r) satisfies the 
quasi-normal boundary conditions, then it can be written 
in the following form 



2Mam 



(r2_ +a2)2^n-i / A'(r+) 
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P(r)=(l-^) V v+ / "^'^^ X J ^ xe^^'^xf(r), (8) 



where F{r) is a regular function at the event horizon 
r = r+ and at the spatial infinity r — oo. 

In order to expand F{r) in the convergent series in the 
region outside the event horizon, we choose the variable 

r — r+ 

z = , 

r - n 

where < r_|_ is chosen so that all the singular points of 
the equation (|4]), except the event horizon and the spatial 
infinity, are located outside the unit circle \z\ > 1. These 
singular points satisfy the equation 

A(r) = 0. 

Fortunately, the structure of the singular points allows us 
to fix for a physically reasonable range of parameters. 
After the value of is fixed, we expand the function 



F into the series 

oo 

F{z)^Y.^,z\ (9) 

Substituting © and © into dU, we find the recurrence 
relation for the series coefficients bi. Following we 
are able to find the equation with an infinite continued 
fraction with respect to w, which can be solved numeri- 
cally. 

The Frobenius method, being based on a convergent 
procedure, allows us to find quasinormal modes (QNMs) 
with any required accuracy, at least theoretically. Prac- 
tically, for some range of parameters (in particularly for 
large azimuthal numbers) the convergence is too slow 
and the calculation cannot be easily done by a personal 
computer. In the range of parameters of slow Frobenius 
convergence, we use the WKB approximation [25^]. The 
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WKB formula ([TO)) was effectively used in a lot of papers 
[29| . and here the WKB accuracy will be checked with 
the help of the the Frobenius method. The 6-th order 
WKB formula reads 

^-gA.=iV+i iV = 0,l,2..., (10) 

where the correction terms A.; were obtained in psj and 
[26j . Here Qo and Qq denote values of Q and its second 
derivative at its maximum with respect to the tortoise 
coordinate r^,, respectively. Since the function fi{uj) is 
found only numerically, in order to find the roots of (|10p , 
we use the following algorithm: 

1. For a given value of w, we find the value of fi using 
the technique described in the section BlI Al 

2. We substitute the numerical values of uj and /i into 
([5]) and find Q as a complex analytical function of 
r. 

3. The analytical continuation of the WKB formula 
into the complex plane implies that "maximum" of 
the potential peak is a complex value, which must 
be correctly chosen among the extrema of the po- 
tential Q. Since the WKB formula provides a suffi- 
cient accuracy only when Re{uj) > Im{uj), we will 
not consider the case, when the imaginary part of 
the quasi-normal frequency is large. Fortunately, 
when the imaginary part of u is small, the imagi- 
nary part of is just a correction to its real part. 
Therefore, the real part of Q{r) is much larger than 
its imaginary part, and the analytical continuation 
of the potential peak is trivial. We choose the loca- 
tion Text of the extremum of Q(r)^ with a real part 
that is larger than the value of the event horizon 
and, at the same time, has the smallest imaginary 
part. This choice of the potential peak does not 
lead to any ambiguities for the considered paramet- 
ric range, is close to the real potential peak for the 
non-rotating case and its searching can be easily 
programmed. 

4. After the potential peak is fixed, we compare the 
frequency iowkb found with the formula ([TU)) with 
our initial guess uj. 

We search for the numerical minimum of the absolute 
value of the difference between the initial guess and the 
WKB formula \ijJwkb{^) — w|, which is a function that 
can be calculated numerically using the above technique. 



IV. QUASINORMAL MODES AND STABILITY 

When we write the quasinormal frequency uj as 

UJ = Re{uj) +ilu\{ijj), (11) 



Im (uj) > corresponds to an unstable (growing) mode, 
and Im (w) < corresponds to a stable (damped) mode. 

The quasinormal modes of a higher-dimensional simply 
rotating black hole are labeled by the quantum numbers 
£, m and j, and also depend on the black hole param- 
eters M and a as well as on the number of space-time 
dimensions D. Thus, the dependence of the spectrum on 
the above six parameters must be carefully considered in 
order to discard or observe instability. We have limited 
the calculations to D — 7,8,9,10 and 11. Keeping the 
event horizon radius r+ fixed, we can measure the QN 
frequencies in units r^^. Thus, the dependence of the 
quasinormal behavior on the four essential parameters 
must be studied. As all of them m, j, a are not lim- 
ited, in principle, we have to check the QN frequency as 
a function of each parameter starting from its minimal 
value and until reaching the asymptotic region where the 
frequency function does not change its behavior qualita- 
tively. However, this task would require us to show here 
hundreds of plots and would take enormous computer 
time. Therefore, we were obliged to adopt the following 
practical strategy. If the frequency a; as a function of 
some parameter does not show any tendency to insta- 
bility and the lowest (longest lived) modes stay very far 
from the instability region, we stopped the calculation at 
some reasonable value of the parameter. Strictly speak- 
ing, this does not allow us to establish mathematically 
rigorous proof of stability, but rather only give a numer- 
ical evidence of stability. In presenting the obtained re- 
sults, in order to reduce the number of plots, we selected 
most representative ones with qualitatively different be- 
havior. If plots are similar (for instance for various D 
and the same other parameters) we have shown only one 
representative plot of the whole "class" . 

Comparison of the accurate data obtained by the 
Frobenius method with the WKB values shows that the 
WKB approximation becomes better as the azimuthal 
number m becomes larger, because it corresponds to 
larger values of the real part of the QN frequency. Indeed, 
the WKB approximation works when |_Re(aj)| 3> |/m(aj)|. 
Thus, the limit of perfect WKB accuracy, which corre- 
sponds to the geometric optics approximation, should be 
achieved at m — >■ cxd. 

Figures [U [2] and |4] show the dependence of the quasi- 
normal frequency, obtained by the Frobenius method, on 
the rotation parameter a for a few low values of m and 
the lowest values of I and j. QNMs for a few low val- 
ues of / and j can be found in Figure [3 There, one can 
see that Re{uj) almost monotonically decreases (except 
for a small peak at very small a), when a is increased, 
while Im(uj) is growing until reaching its maximum, and 
then falling off to some asymptotic constant value. As 
the absolute value of Im{uj) decreases when m increases, 
it is important to find QNMs for sufficiently large m, in 
order to rule out the instability. On figures Fig. [6] one 
can see the high m behavior for a fixed a and D. Ap- 
parently Re(uj) is equidistant in the limit m — > oo. The 
absolute value of Im{uj) approaches some maximum and 
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FIG. 2: Quasinormal modes obtained by the Frobenius method for D = 7, I = 2, j — 0: m — —1 (blue), m — —2, m = —3, 
m = —4, m = —5 (hght blue), m = —10 (green), m — —15 (red), m = —20 (magenta), m = —25 (brown). Higher negative 
values of m correspond to the larger real and imaginary part of to. The imaginary part of the QN frequency stays negative, 
implying stability against perturbations with the high negative azimuthal number m. 
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FIG. 3: Maximum value of the imaginary part of the quasi-normal mode for negative values of m together with the found fit 
max{Im{uj) cc (— m)~" (a > 0) for D = 7,l — 2, j — (left figure, log-log scale) and corresponding values of the rotation 
parameter Umax as a function of — m (right figure). 



6 




FIG. 5: Quasinormal modes obtained by the Frobenius method for D = 7: {I — 2, m — 0, j — - black), {I — 2, m — 0, j — 1 

- yellow), {I = 2, m — 0, j = 2 - orange), (Z = 2, m = 1, j = - purple), (l = 2, m = 1, j = 1 - blue), (Z = 2, m = 1, j = 2 

- magenta), {I — S, m = 1, j — - green), {I = 4, m — 1, j = - cyan), {I — 5, m = 1, j = - red). The larger values of j 
correspond to the smaller values of the imaginary part of to. The real part increases with l + m + 2j for small rotation. When a 
is large, the real part for the multipole index j ^ approaches zero, while that for j — approaches a constant that increases 
with I. 
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FIG. 7: Quasinormal modes as a function of a for D — 8, I — 2, m = 0,1: j = (red), j — 1 (green), j = 2 (blue). Higher m 
and j have larger real and imaginary part of cj 



Re{w) 




FIG. 8: Quasinormal modes as a function of a for D — 9, I — 2, m = 0,1: j = (red), j — 1 (green), j = 2 (blue). Higher m 
and j have larger real and imaginary part of to. 



then starts to fall off at large m, reaching some constant 
value asymptotically (see Fig. |6|). 

In the figure[2]one can see that for large negative m the 
maximum of the imaginary part is higher (closer to zero) 
and occurs for larger value of the rotation parameter a. 
One might ask himself if there is sufficiently high value 
of m for which the imaginary part crosses zero at some 
point leading to the instability for larger values of a. Our 
calculations show that this does not take place. In the left 
figure 13] we plot the dependence of the maximum value of 
Im (uj) as a function of —m in the logarithmic-logarithmic 
scale. We see that this value is well approximated by the 
potential law maa;(Im (w)) cx (— m)"" (a > 0). If this is 
valid for the asymptotically large values of m, the maxi- 
mum of the imaginary part of the quasi-normal frequency 
approaches zero only asymptotically. The corresponding 
value of the rotation parameter depends linearly on m 
and, therefore, also approaches infinity (see the right fig- 
ure [3|). Thus, we conclude that at least for the finite 
rotation parameter all modes of negative m are stable. 

In figures H [51 d [TO] the dependence of QNMs on j 
is shown. There, one can see that the QN spectrum 
stays well far from the instability point. Thus, we con- 



clude that no growing mode was found for simply rotating 
Myers- Perry black holes with D = 7,8, 9, 10 and 11, i.e. 
tensor-type gravitational perturbations are stable in the 
above cases. 



V. CONCLUSIONS 

We have investigated the stability of simply rotating 
Myers-Perry black holes with D > 7 against tensor-type 
gravitational perturbations. We have not found a grow- 
ing mode, which apparently signifies stability against ten- 
sor perturbations. 

In addition, we gave here extensive numerical data for 
the quasinormal frequencies for the tensor-type gravita- 
tional perturbations of these black holes. 

Let us note that the quasinormal spectrum of a simply 
rotating Myers-Perry black hole is not a limit of that of 
a MP black hole in the AdS spacetime, when the AdS 
radius R approaches infinity. From the equation (71) of 
[2l| . one can see, that when R ^ r+ the quasinormal 
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Re(w) 




FIG. 9: Quasinormal modes as a function of a for D — 11, I — 2, m = 0,1: j — (red), j = 1 (green), j = 2 (blue). Higher m 
and j have larger real and imaginary part of uj. 
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FIG. 10: Dependence of the quasinormal modes on j for D = 8 (blue) and D = 9 (red), a — 7, m = 1. The plots with larger 
Re{uj) and Imiio) correspond to the larger D. 



modes approach zero as 

2n + D + 2l + m + j- l + F{D,l,m,j, a) 



R 



= O 



In the forthcoming paper [3011 , we shall consider gravi- 
ton emission in the bulk for simply rotating black holes 
as well for black holes with all equal angular momenta. 



ginning of this work R. A. K. was supported by the Japan 
Society for the Promotion of Science (Japan) and at the 
final stage by the Alexander von Humboldt Foundation 
(Germany) R. A. K. also acknowledges the hospitality of 
the Theory Division of the High Energy Accelerator Re- 
search Organization (KEK) in Tsukuba, where a part of 
this work was done. A. Z. was supported by Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), 
Brazil. 



Acknowledgments 

H. K. was supported in part by Grants-in-Aid for Sci- 
entific Research from JSPS (No. 18540265). At the be- 



[1] R. Emparan and H. S. Reall, Living Rev. Rel. 11, 6 [3] H. Kodama, Lect. Notes Phys. 769, 427 (2009) 

(2008). ;arXiv:0712.2703 [hep-th]]. 

[2] R Kanti. larXiv:0802. 22181 [hep-th]. [4] Kodama, H.: Prog. Theor. Phys. Suppl. 172, 11-20 



9 



(2008). 

[5] Kodama, H., Ishibashi, A. and Seto, O.: Phys. Rev. D 

62, 064022 (2000). 
[6] H. Kodama and A. Ishibashi, Prog. Theor. Phys. Ill, 

29 (2004). 

[7] A. Ishibashi and H. Kodama, Prog. Theor. Phys. 110, 
901 (2003). 

[8] H. Kodama and A. Ishibashi, Prog. Theor. Phys. 110 
(2003) 701. 

[9] R. A. Kono plya and_A. Zhidenko, Nucl. Phys. B 777, 
182 (2007) arXiv:hel>th/0703231 . 

[10] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 78, 
104017 (2008) arXiv:0809.2048 [hep-th]]. 

[11] R. A. Konoplya and A. Zhidenko, ar Xiv:0809.2822 [hep- 
th]. 

[12] T. Takahashi and J. Soda, larXiv:09 02. 2921' [gr-gc]. 

[13] M. Beroiz, G. Dotti a nd R,. j. Gleiser, Phys. Rev. D 
76, 024012 (2007) arXiv:iiep-th70703074 ; R. J. Gleiser 
and G. Dotti, Phys. Rev. D 72, 124002 (2005) 
(arXiv :gr-qc/0510069 . 

[14] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 77, 
104004 (2008) arXiv:080 2.0267 [hep-th]]. 

[15] M. Kimura, K. Murata, H. Ishihara and J. Soda, Phys. 
Rev. D 77, 064015 (2008) arXiv:0712.4202 [hep-th]]; 
H. Ishihara, M. Kimura, R. A. Konoplya, K. Murata, 
J. Soda and A. Zhidenko, Phys. Rev. D 77, 084019 (2008) 
[arXiv :0802.0655 [hep-th]]. 

[16] R. A. Konoplya, K. Murata, J. Soda and A. Zhidenko, 
Phys. Rev. D 78, 084012 (2008) arXiv:0807.1897 [hep- 
th]]. 

[17] S. R. Dolan, Phys. Rev. D 76, 084001 (2007) 

[arXiv :0705.2880 [gr-qc]]. 
[18] H. K. Kunduri, J. Lucietti and H. S. ReaU, Phys. Rev. D 

74, 084021 (2006) arXiv:hep-th/0606076 . 
[19] K. Murata and J. Soda, arXiv:080 3.1371. [hep-th]. 
[20] K. Murata, arX iv:0812.0 718 [hep-th]. 
[21] H. Kodama, R. A. Konoplya and A. Zhidenko, 



arXiv:0812.0445 [hep-th]. 

[22] H. Suzuki, E. Takasugi a nd H. U metsu, Prog. Theor. 
Phys. 100 (1998) 491 arXi v:gr-qc/9 805064 ; R. A. Kono- 
plya and A. Zhidenko, Phys. Rev. D 76, 084018 (2007) 
arXiv:0707.1890 [hep-th]]. 

[23] E. W. Leaver, Proc. Roy. Soc. Lond. A402, 285 (1985). 

[24] A. Zhidenko, Phys. Rev. D74, 064017 (2006), 
gr-qc/0607133 ; R. A. Konoplya and A. Zhidenko, Phys. 
Rev. D 73, 124040 (2006) arXiv:gr-qc/0605013 . 

[25] B. F. Schutz and C. M. WiU Astrophys. J. Lett 291 L33 
(1985); S. Iyer and C. M. WiU Phys. Rev. D 35 3621 
(1987). 

[26] R. A. Konoplya, Phys. Rev D 68, 024018 (2003); 

R. A. Konoplya, J. Phys. Stud. 8, 93 (2004). 
[27] R. A. Konoplya, Phys. Rev. D 68, 124017 (2003) 



[28] 



arXiv:hep-th/0309030t . 



H. T. Cho, A. S. CorneU, J. Doukas and W. Naylor, 
arXiv:0904.1867 [gr-qc] . 
[29] E. AbdaUa, O. P. F. Piedra and J. de Oliveira, 
larXi7:0 810.5489 [hep-th]; M. 1. Liu, H. y. Liu and 
Y. X. Gui, Class. Quant. Grav. 25, 105001 (2008) 
[arXiv: 0806.2716 [gr-qc]]; P. Ka nti and R. A . Konoplya, 
Phys. Rev. D 73, 044002 (2006) 'arXiv:hep-th/0512257; ; 
J. F. Chang, J. Huang and Y. G. Shen, Int. J. Theor. 
Phys. 46, 2617 (2007); H. T. Cho, A. S. Cornell, 
J. Dou kas and W. Na ylor, Phys. Rev. D 77, 016004 
(2008) [arXiv:0709. 16611 [hep-th]]; J. F. Chang and 
Y. G. Shen, Int. J. Theor. Phys. 46, 1570 (2007); 
Y. Zhang and Y. X. Gui, Cla ss. Quant. Grav. 23, 6141 
(2006) arXiv:gr-qc/0612009"; R. A. Konoplya, Phys. 
Lett. B 550, 117 (2002) arXiv:gr-qc/0210105 ; S. Fer- 
nando and K. Arnold, Gen. Rel. Grav. 36, 1805 (2004) 
arXiv:hep-th/031204l]. 

P. Kanti, H. Kodama, R. A. Konoplya, N. Pappas and 
A. Zhidenko. .arXiv:0906.3845l [hep-th] . 



[30] 



